Magnetization plateau and quantum phase transitions in a spin-orbital model 
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A spin-orbital chain with different Lande g factors and one-ion anisotropy is studied in the 
context of the thermodynamical Bethe ansatz. It is found that there exists a magnetization plateau 
resulting from the different Lande g factors. Detailed phase diagram in the presence of an external 
magnetic field is presented both numerically and analytically. For some values of the anisotropy, the 
four-component system undergoes five consecutive quantum phase transitions when the magnetic 
field varies. We also study the magnetization in various cases, especially its behaviors in the vicinity 
of the critical points. For the SU(4) spin-orbital model, explicit analytical expressions for the critical 
fields are derived, with excellent accuracy compared with numerics. 



Orbital degeneracy in electron systems leads to rich 
and novel magnetic phenomena in many transitional 
metal oxides [1]. Among them are the orbital ordering 
and orbital density wave, which have been observed ex- 
perimentally in a family of manganites [2]. A tractable 
model to describe 2-fold orbital degenerate system is the 
SU(4) model [3], which has attracted much attention 
[3-9]. In the one-dimensional case the model is exactly 
solvable by Bethe ansatz (BA) [4]. An interesting ques- 
tion is to study the critical behavior of such a system 
in an external magnetic field, especially when different 
Lande g factors for spin and orbital sectors are involved. 
One may expect that the difference of g factors will bring 
about new physics as a result of the competition of the 
spin and orbital degrees of freedom. In Ref. [5] , the au- 
thors studied the magnetic properties of the SU(4) model 
via BA, without taking different g factors into account, 
whereas numerical calculation was performed in Ref. [6] 
for the model with different g factors for up to 200 lat- 
tice sites. However, a full picture about the critical fields 
is still lacking. Another motivation is to see whether or 
not any magnetization plateau (MP), an interesting mag- 
netic phenomenon, occurs in such a spin-orbital model. 
As is well known, antiferromagnetic chains with integer 
spin are gapful [10], whereas for half-integer spin there 
also exists a gapful phase with a MP in the presence 
of a large planar anisotropy [11]. Also fractional MP 
have been observed and can be explained by Shastry- 
Sutherland lattice [12]. But an MP arising from different 
Lande g factors has not been addressed yet. 

Deviation from the SU(4) symmetry was considered by 
variation in the interaction parameters of neighbor sites 
[7,8], while another possible deviation may result from 
the one-ion interaction. Since many compounds are mag- 
netically anisotropic in which the orbital angular momen- 
tum (0AM) may be constrained in some direction due to 
crystalline field, the angle between spin and 0AM deter- 
mines the spin-orbital coupling (SOC) energy. This kind 



of one-ion SOC leads to magnetic anisotropy [13]. Un- 
der the influence of molecular field and an external field, 
the spin is parallel to the 0AM, then the SOC energy de- 
pends on whether they are in the same directions. In such 
a case, sfrf type of interaction describes well the SOC 
energy. Here we introduce such an SU(2)®SU(2) SOC 
interaction into the SU(4) model. A detailed investiga- 
tion of the phase diagram is undertaken both numerically 
and analytically in the context of the thermodynamical 
Bethe ansatz (TBA). We find that the system exhibits 
an MP resulting from different g factors when the SOC 
is sufficiently strong. The critical behavior of the magne- 
tization in the vicinities of the critical points is revealed. 
For certain values of Lande g factor, the model undergoes 
five consecutive quantum phase transitions when the ex- 
ternal magnetic field varies. Further, the explicit analytic 
expressions for the critical fields for the SU(4) model are 
derived, with excellent accuracy compared to numerical 
results. 

The model and TBA. Consider an L-site chain with the 
Hamiltonian 

1=1 

Hz = AzY^-str,^, M ^ ~g,Hj2st ~ gtHj^Tl, (1) 

i i i 

where s and r are spin- 1/2 operators for spin and 
orbital sectors, and gs and gt denote, respectively, 
the Lande g factors. We assume gs > gt through- 
out the paper. Ho is the SU(4) model with Pij — 
{2si ■ Sj -\- 1/2) {2fi ■ fj + 1/2) exchanging the four site 
states [sfrf ): 0i = [ti), 02 = lit) , 03 = ITT), 04 = |ii). 
The symmetry is broken into SU(2)(g)SU(2) by Hz and 
further into four U(l)'s by the external magnetic field 
H. The model can be solved exactly via BA approach. 
The BA equations are the same as the SU(4) model [4,14] 
under the periodic boundary conditions, with the energy 
eigenvalues given by 
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where a„ (A) = l/(27r) n/(A2 + nV4), and Ei = 
-A,/4 - 5_i//2, ^2 = -A,/4 + g_H/2, E3 = A,/4 - 
g+H/2, Ei = + g+H/2, with 5(± = 5^ ± g*. A^fe is 
the total site number in state ^fc and M^'^ [i = 1,2,3) 
is the rapidity number. For a certain choice of the basis 
order, which depends on whether or not the component 
is energetically favorable, the energy can be rewritten 

as E = Y.ti g^^H^i) + g^^^M^^) + ff(3)M(3), with the 
coefficients g*^*' depending on the basis order. In the ther- 
modynamical limit, the summation over the rapidity be- 
comes integral. Following Refs. [15] and [16], one may 
obtain the ground state (GS) equations for the dressed 
energies e^'' {i — 1,2,3), 
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where e^°^ = e^^^ = and the symbol * denotes the 
convolution. The dressed energies describe elementary 
excitations over the Fermi seas, i.e., the GS with all 
states with negative dressed energies filled. According 
to an energetics argument, we may divide the exter- 
nal field H into three regions: (I) < < Hrx, 
(II) Hri < H < Hr2, (III) Hr2 < H < 00 with 
Hri = |A,|/(2g,), Hn2 = I A, 1/(2.9,). For A, > 0, 
the corresponding basis order arc: (1+) (01,(562, 0;^, 04)"' , 
(11+) i^,,^3,h,^4V, (HI) ih,^i,4>2,4>4V; for A, < 0: 
(I_) {h,'p4,<l^i,h )^, (II_) {h,4'i,^4,hV, (III) the 
same as A^ > 0. These five basis orders provide a full 
description of the phase diagram of the system. 

Plateau. The competition between the anisotropy pa- 
rameter Az and the magnetic field H results in a novel 
quantum phase diagram. In the absence of the magnetic 
field, it is easy to find that the states 03 and 04 are gap- 
ful for A^ > A^ = 4 In 2. Whereas for A^ < -A^, the 
components 0i and 02 are gapful. Therefore, the GS is 
in an su{2) spin-orbital liquid state in strong anisotropy 
regime in the absence of the field. However, the presence 
of the magnetic field completely splits all foiu compo- 
nents energetically. The magnetization — ggS^ + gtT^ 
increases from zero. For large positive A^, the field bring 
the component (j)3 closer to the GS, while the component 
4>2 gradually gets out of the GS. Certainly, if the field 
reaches the first critical field where the component (^3 
has not yet involved in the GS, a quantum phase tran- 
sition from the spin-orbital liquid phase to a ferromag- 
netic phase occurs. Thus a magnetization plateau opens 
with a constant magnetization = (7-/2. Neverthe- 
less, this plateau will disappear when the field is strong 
enough H > H^2- such a case, the component (^3 
becomes involved in the GS. The critical field H^2 indi- 
cates a quantum phase transition from the ferromagnetic 
GS into a spin-orbital liquid phase. If the field continues 



to increase, the spin and orbital sectors become fully- 
polarized at the third critical point H^^. From the TBA 
equations (3), we get the exact expressions for the critical 
fields 



rp — — rp 



A,/2 - 4 _ A,/2 -I- 4 



(4) 



Notice that the plateau opens only if A^ > A^ = 8gs/g- 
and < gt < gs- If the g factors are the same, the 
plateau disappears because the components 0i and 02 
remain degenerate in the field. The critical behavior of 
the magnetization in the vicinities of the critical points 
may be summarized as follows 



(5) 



where (M^)^ = gz^/2 are, respectively, plateau and satu- 
ration magnetizations. 5H is the small deviation from 
the critical points and ry = ±1 depending on is 
increasing or decreasing. Um = g^^/T^ near H^-^ and 

km = gV^ /"^ near H^^ ^-^id ^c3- "^he coefficients g^ and 
gt in (4) and (5) can be easily understood, since only 
and 02 exist in the GS before H^^ is reached, the dif- 
ferences of their field energy and magnetization are g-H 
and g- respectively. While for i/^z < H < H^^, only 
(j)! and compete in the GS, since 02 already gets out 
before (^3 enters the GS. (l>i and 03 differ in field energy 
by gtH and in magnetization by gt- In Fig.lA we plot 
the magnetization curves for different values of the pa- 
rameter A2, including a plateau case. 

The phase diagram: numerical and analytical. Here 
we present a detailed analysis of the GS phase diagram 
both numerically and analytically. In Fig. 2, we plot the 
phase diagram with respect to the parameters A^ and 
H for fixed values of gs and gt {gs = 2.0 and gt = 1.0). 
For convenience, we refer to the GS with i components 
as i-state GS. Then for the phase transition between 3- 
state and 2-state GS, the critical fields follow from the 
Wiener-Hopf method [17], which is valid for large Fermi 
boundaries. Explicitly, we have 

H^^ = (A, - A^)5;i - T,g'_g^'{Az - A^)^ 
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29s - gt 



[ff«A g-(ff, + g,)A.]^ 

(2y.s+y/)-'' 



A'' - ^ 
^gs + gt 

H^^- = 7r'g-gf[^l - nAgZ'gliA, + Ao)-l], (6) 



where n = l/(27r2) and Af = 2gtH^ = gt/g-Al which 
corresponds to infinite Fermi boundaries. Point Q is 
determined by Aj = -2gtH'^ , with H'^ = 2\a2/gs + 
47r~^ In^ 2gt' lg\. Near P or F the above analytic results 
deviate due to small Fermi boundaries. But in this case, 
an analysis may be carried out in terms of expansion of 
small Fermi boundaries. This leads us to H^^ , Hc^* = 



2 



1^ ± (5-/2)i/(7rff|)(Af - A,)i. Similarly, near F, we 

have HfQ, = -I: ± (.9+/2)i/(^,9l)(A, - Af )i. 

For lA^I < A°, the GS involves all the four com- 
ponents, the magnetic field first brings about a phase 
transition from a 4-component liquid to a 3-component 
liquid at the phase boundary C+QC_. Here one of 
the four components, which is energetically unfavorable, 
completely gets out of the GS, then the corresponding 
Fermi sea disappears. So the critical field only involves 
two Fermi points Bi and B2. At point C+, the compo- 
nent <j)2 is degenerate with the component (pi which is 
the most energetically favorable, the Fermi point Bi lies 
at infinity. Increasing H along C+MNQ drives (t)2 away 
from so (j)2 becomes less energetically favorable in 
the GS. Therefore the first Fermi sea shrinks, i.e., Bi de- 
creases from infinity. Beyond M point, both the increase 
of H and the decrease of A^ make; (/>5 sink below the (j)2 
which is rising, the energy difference between (j)^ and (/>i 
begins to dominate over Bi. As ^3 is drawing near ^1, 
the first Fermi sea becomes broadened again with an in- 
crease of Bi. After point N, ^3 sinks beyond (/>i to be 
the lowest state, <?!>i becomes less favorable in the GS. 
The first Fermi sea shrinks again, Bi begins to decrease 
along NQ from the infinity at N. A similar analysis is ap- 
plicable to B2 , which is mainly influenced by the energy 
difference between the second and third components in 
energy levels. B2 rises from zero at C+ to infinity at M 
and decreases along MNQ to zero at Q. In the respective 
sections of C+MNQ the critical fields take the form 

jjC^MNQ ^ ^^^^ig^ ^ jjC+MNQ ^ ^jjC^MNQ 

H+ + T2/^l,Jg+ + wT2{g-ao - .9tA„,3)V3+ 
^= H+ + T2[{g-ao + g.A„,3)2 + w{g-aa - gtAa,3f]/gl 
^= H+ + T2{g-ao + (7,A„,3)Vff+ + wr2AlJg+, 

where Aa,m = o-o - ao = ^ti" - |ln3, r2 = 

and w — 2/3. In each case, the first term comes from 
infinite Bi and B2, the second term is correction from 
finite but large B2 , and the third term is the leading cor- 
rection from the larger Bi. For C+MN near point M 
the B2 and Bi terms need to be exchanged since B2 be- 
comes larger. Along NQ Bi is always larger than i?2- 
The location of Bi = B2 in C+MN may be estimated 
by Hb^ = Hb2, which gives A^ = 2aa{2gs ~ gt)/{igs) for 
C+M and A^ = |ao for MN. This coincides well with nu- 
merics, e.g., for MC+ and gg =2.0, gt =1.0, the analytic 
result is A^ =1.073 whereas the numerical one is 1.042. 

Similarly, for C_VQ, the term resulting from the in- 
finite Fermi boundaries is = Aa^-i/g-, the cor- 
rection terms are H^^ ~ = T2A^ _^/g_, H^^ ~ = 

Hb? = Mg+ao + gtAa,-3?/gt and H^^ = rs^+ao + 
5sAa^_3]^/gi, respectively. These expressions are not 
valid for gs ~ gt due to small Fermi boundaries. 



When the system is fully-polarized, only (f)^, exists in 
the GS, while the other components are all gapful, with 
a gap A = mm{Ei - 4 — ^^3 | i = 1,2,4}. This gap is 
closed if i? < Hf, with the fully-polarized critical point 

Hf = max{{A,/2 + 4) g^\ Ag-'}. (7) 

This expression is exact and valid for all A2. When 

Az < Af = —8gs/g+, the strong negative anisotropy 
makes and (j)2 too far away from (j)^. Before the field 
brings them close enough to get involved in the GS, the 
component (f>4 has been all pumped out by the field from 
the GS at critical point Hf = 4:/g+. For all A^ < Af , 
the magnetization is the same as shown by curve a in 
Fig.lA. 

The analytic results are compared with the numerics 
in Fig. 2, with very satisfactory accuracy. 

Five consecutive He phase transitions. The competi- 
tion of anisotropy A^ and the external field H also leads 
to an unusual magnetic phenomenon. For some fixed val- 
ues of Az, gs and gt, the system undergoes five consecu- 
tive quantum phase transitions when H varies, although 
the 4-component model usually has at most three consec- 
utive phase transitions. A strong negative anisotropy A^ 
makes 04 energetically quite favorable. The field H ex- 
pels (1)2 first from the 4-component GS. However, before 
H overwhelms the influence of A^ on (^4, further increase 
of H will make 4>2 closer to and draw it back into 
the GS. This process brings about the first two phase 
transitions. Then H plays a dominant role, it begins 
to bring out (/>4, 4*2 and (/>i from the GS one by one. 
This results in other three consecutive phase transitions. 
The variation of the state component numbers in the GS 
is: 4^3^4^3^2^1. This five H^s case exists for all 
< gt < gs and becomes more visible when gt is larger. 
One case is marked by black dots in Fig.3B for gs = 2.0 
and gt = 1.8. Another possible case for five He transi- 
tions to occur is the GS composed of ^3 and 04. The 
field brings (pi into the GS first. As the difference of 
gt and (/« is getting smaller, 02 has closer energy to 0i. 
Further increase of H brings 02 into the GS before 04 
completely gets out. The component number changes in 
such a way: 2^3^4^3^2^1. This case occurs when 
the point Q is below C_ in Fig. 2, approximately requir- 
ing 1 > gt/gs > 1 — 2ln2gf /{n'^g'^). Four consecutive 
He transitions take place for A^ = — A^. In Fig. IB and 
C, we plot the magnetization curves which display five 
consecutive He phase transitions. 

The SU(4) model. If we set the anisotropy parameter 
A^ to be zero, the model reduces to the SU{4:) model 
with different Lande g factors in the spin and orbital sec- 
tors. In this special case, the above results for the critical 
fields give rise to 

Tj-SU(4) ■ -I , 2/0 \2 -3 , 2 -1 

Hei =aog_^ +T2aQ(2gs - gt) g+ +wT2aQg_^_ , 
jjSUiA) ^ ^e ^2^^ _ ^^^-1 ^ ^^^g^Ain2gs - gt)-\ (8) 
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We compare the above analytic results with TBA nu- 
merical ones in the Fig. 3A, which shows an excellent 
accuracy for most values of gt- Notice that the deviations 
occur only for less-physical values of gt due to small Fermi 
points. For gg = 2.0 and gt = 1.0, the discrepancies of 
H^^^^^ and H^^'^^^ (TBA numerics: 0.3695, 0.9415; ana- 
lytic: 0.3697, 0.9386 (200 sites results [6] for comparison: 
0.31, 0.93)) are, respectively, 0.05% and 0.3%. ^£3'^^"^ is 
exact. The magnetization of the SU(4) model is shown 
in Fig.lA, which also coincides with the numerical result 
for 200 sites [6]. 
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FIG. 1. (A) Typical magnetization behaviors for fixed values of Lande g factors Qs = 2.0 and gt = 1.0. The numbers in 
brackets indicates values of the anisotropy parameter A^. The dotted line e denotes the magnetization for the SU(4) model. 
Curve j exhibits a magnitization plateau. (B) Magnetization corresponding to five consecutive five He quantum phase transitions 
for gs — 2.0, gt ~ 1.9, and Az = —3.0. Here the variation of the state component numbers is 2^3 — >4— >3^2— >1. The number 
i labels the state 0i, e.g., the state components in the phase 123 are 4'i4>24>i in which 4>i is the most energetically favorable 
whereas 03 is the least energetically favorable. The phase variations between Hc2 and Hc3 are 3412^3142^3124. The dotted 
line for comparison is an extension of the phase 34 by assuming the components unchanged. (C) Magnetization for consecuitive 
five He phase transitions for gs = 2.0 and gt = 1.0, with Az = —1.4. Hci =0.531 (transition 3142—^214) and Hc5 ~ 3.33 
(31^3) are relatively far away. The variation of the component numbers in the phase transitions is 4-^3^4^3^2^1. 
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FIG. 2. Phase diagram for fixed values of Lande g factors ga — 2.0 and gt ~ 1.0. Inside C+QC_, the 4-state phases are 1234 
(I+), 1324 (II+), 3124 (III), 3142 (11-), and 3412 (I_). The discrepancy of the analytic curves from the numerical ones is not 
visible for most regions of C+PQ and C+MNQ (with typical differences within 1.0%, and 0.1%, respectively). There is less 
accuracy in regions I_ and II_ for the analytic results due to larger driving terms and smaller Fermi boundaries. Magnetization 
plateaux and fully-polarized cases are exact. 
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FIG. 3. (A) Comparison of the analytic and the numerical results for the three critical fields of the SU(4) model with respect 
to the orbital g factor gt for a fixed value of the spin g factor gs — 2.0. The numbers label the state components as in Fig. 2, 
and Hc3 (phase transition 31— >3) is exact. (B) Phase diagrams for various gt {gs = 2.0). = corresponds to the SU(4) 
model. The black dots mark consecutive five He phase transitions. Such a phenonemon exists for all < gt < gs, with Hd-Hd 
being closer for smaller gt- 
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